Possible Evidence For Axino Dark Matter In The Galactic Bulge 



Dan Hooper^ and Lian-Tao Wang^ 
^Astrophysics Department, University of Oxford, Oxford, UK; ^Department of Physics, University of Wisconsin, Madison, 

USA. 
(February 7, 2008) 



^ ■ 

o ■ 
o : 

(N . 

> : 
o . 

o\ : 

m ■ 
> 
o 

(N 
(N 
(N 
O 

o : 

i!^ ■ 
Oh; 

I 

Oh! 



X 



Recently, the SPI spectrometer on the INTEGRAL sateUite observed strong 511 keV line emission 
from the galactic bulge. Although the angular distribution (spherically symmetric with width of 
~ 9°) of this emission is difficult to account for with traditional astrophysical scenarios, light dark 
matter particles could account for the observation. In this letter, we consider the possibility that 
decaying axinos in an R-parity violating model of supersymmetry may be the source of this emission. 
We find that ~ 1 — 300 MeV axinos with R-parity violating couplings can naturally produce the 
observed emission. 
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I. INTRODUCTION 

The SPI spectrometer on the INTEGRAL (INTEr- 
national Gamma-Ray Astrophysics Laboratory) satel- 
Ute has made the observation of a bright (9.9^21 x 
10~^ phcm~^ s~^), 511 keV gamma-ray emission line 
from the galactic bulge [1]. The emission is consistent 
with being spherically symmetric, with a full- width-half- 
maximum of about 9° (6° — 18° at 2a confidence). The 
3 keV width of the line is dominated by e -I- e~ annihi- 
lations via positronium formation [2] . These findings are 
in agreement with earlier observations, such as those by 
the OSSE experiment [3]. 

This observation of bright 511 keV emission from the 
galactic bulge has been quite difficult to explain with tra- 
ditional astrophysics. Most potential sources considered 
do not produce a sufficient number of positrons (such 
as neutron stars, black holes [4], radioactive nuclei from 
supernovae, novae, red giants or Wolf-Rayet stars [5], cos- 
mic ray interactions with the interstellar medium [6] , pul- 
sars [7] or stellar flares) and those which may possibly be 
capable of producing the required number (type la su- 
pernovae [8,9] or hypernovae [10,9]), may not be capable 
of filling the entire galactic bulge [11]. In particular, the 
rate at high altitude is likely to be too low to explain the 
observed extension of the 511 keV source [12]. 

Given this difficulty, alternative explanations should 
be considered. In Ref. [13], it was suggested that light 
(1-100 MeV) dark matter particles annihilating in the 
galactic bulge could explain the observed emission. In 
this letter we, instead, consider light decaying dark mat- 
ter particles. In particular, we consider the supersym- 
metric partner of the axion [14], the axino, in R-parity 
violating supersymmetric models. 

In the Minimal Supersymmetric Standard Model 
(MSSM), a Z2 symmetry, R-parity [15], is usually im- 
posed to forbid dimension four operators which lead to 
fast proton decay [16]. A by-product of an exact R-parity 
is that the Lightest Supersymmetric Particle (LSP) is 
stable and often a good candidate for cold dark mat- 



ter. Although R-parity is an elegant way of suppressing 
proton decay, it does not have to be the only way that 
nature can choose to do so. In particular, baryon parity 
is sufficient to make the proton stable. In this case, the 
coupling strengths of other R-parity violating operators, 
such as XijkLiLjE^, are much less constrained. For a 
summary of R-parity violating coupling constraints, see 
Ref. [17]. 

A typical dark matter candidate provided by R-parity 
conserving supersymmetry is a neutralino LSP. Neutrali- 
nos are the superpartners of the neutral gauge bosons 
and Higgs bosons and have masses constrained by di- 
rect searches to be greater than ~ 30 GeV, too heavy 
to produce a large fiux of thermal positrons. Of course, 
with sizable R-parity violation, the neutralino will have 
a short lifetime and cease to be a good candidate for cold 
dark matter. 

With the Peccei-Quinn (PQ) mechanism [14] providing 
a natural solution to the strong-CP problem in a super- 
symmetric theory, the existence of an axino is inevitable. 
The axino's mass is expected to be considerably lower 
than the electroweak scale, perhaps in the keV to several 
GeV range and is capable of providing the observed quan- 
tity of dark matter given a low reheating temperature 
[18-22]. Since the axino is in the same supermultiplet as 
the axion, its couplings to matter fields are generically 
suppressed by /~^, where fa is the PQ symmetry break- 
ing scale ~ 10® — 10^^ GeV. As we shall see in detail in 
this paper, due to this large suppression, axinos could 
be considered stable during the history of the universe 
even in the presence of sizable R-parity violating cou- 
plings. Furthermore, a long-lived MeV-GeV axino, such 
as we consider in this letter, would be heavy enough to 
constitute a good candidate for cold dark matter [21-23]. 

A possible concern for axinos in the early universe 
is their effect on the light element abundances [24-26]. 
With a long-lived axino, however, axino decays do not 
threaten these observations. Furthermore, with the in- 
troduction of R-parity violation in our scenario, we al- 
low for relatively fast Next-to-Lightest Supersymmetric 
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Particle (NLSP) decays into Standard Model particles. 
Thus the epoch of SUSY particle decays is over prior to 
nucleosynthesis and the light element abundances remain 
unaffected [27]. 

Before going into a more detailed study of decay- 
ing axinos, we briefly consider here the possibility of 
decaying gravitinos. Gravitinos also have highly sup- 
pressed couplings to matter {Mp^ in the kinemati- 
cal regime of interest) and could be a good candidate 
for cold dark matter. We have found, however, that 
for the case of trilinear R-parity violating terms, the 
gravitino lifetime is too long to account for the ob- 
served 511 keV emission (see section II). In particu- 
lar, the gravitino lifetime is estimated to be T3/2 ~ 
10^1 sec (mj/lOOGeV)'' (O.l GeV/ma/a)^ {0.1/Xf. 



II. DECAYING DARK MATTER 

If light dark matter particles constitute the galac- 
tic halo, decays of such particles can produce positrons 
which eventually annihilate producing the observed 511 
keV emission. In this section, we calculate the lifetime of 
a light decaying dark matter particle needed to account 
for the observed 511 kcV flux. 

We note that although the observed emission line is 
energetically narrow, indicative of positrons annihilating 
at rest, positrons resulting from dark matter decays need 
not be at rest, as they will stop easily via ionization losses 
before annihilating as long as their initial energy is less 
than ^ 100 MeV [13]. This energy corresponds to de- 
caying dark matter particles of a few hundred MeV or 
less. 

The dark matter distribution in the galactic halo is 
traditionally parameterized by 
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{r/ap[l + ir/ap](0-^)/»' 



(1) 



where a, (3 and 7 are given by the choice of halo profile 
and o is the distance from the galactic center at which 
the inner power law breaks. In the galactic bulge, r -C 
a, and the parameterization reduces to p{r) oc (j-Ja)~i ■ 
Integrating over the line-of-sight of the observation, and 
averaging over the angular resolution of SPI 2°), the 
angular distribution of 511 keV gamma-rays from dark 
matter decays for a given halo profile can be compared to 
the observations of SPI/INTEGRAL. In figure 1, we show 
these results. We find that for a cusped profile, 7 ~ 1.2, 
the data is well fitted. This is in agreement with the 
results of recent high resolution N-body simulations [28] . 
Note that in Ref. [13], the best fit was found for 7 ~ 
0.6. This is because for annihilating dark matter, the 
annihilation rate is proportional to the density squared, 
rather than simply the density. 
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FIG. 1. The angular distribution of 511 keV 7-rays from 
decaying dark matter averaged over the 2° angular resohition 
of the SPI spectrometer on INTEGRAL for several halo pro- 
files. SPI's observation indicates a full width half maximum 
of 9° with a 6° — 18° 2cr confidence interval. Shown as vertical 
dashed and dotted lines are the central value and 2cr limits of 
the angular widths found by SPI. To agree with this data, a 
cusped halo model with 7 ~ 0.8-1.5 is favored. 

Normalizing the halo profile to the local dark matter 
density, we obtain 



p{r) ~ 



0.3Mq/pc3 
(r/lkpc)'^ ' 



(2) 



where 7 ~ 1.2. The total mass within the 9° circle ob- 
served by INTEGRAL is then 
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670 pc 



p{r)A-Kr^dr ~ 1 x 10^ Mg ^ 1.5 x Iff" GeV 



(3) 



The rate at which dark matter decays in this region is 
simply the total mass divided by the mass and lifetime 
of the dark matter particle. Also, the rate of decays can 
be matched to the observed flux of 511 keV gamma-rays. 
Comparing these two quantities yields 



1.5 X 1066 GeV 

^dm'^dm 



7,51l47ri?Qc- 



(4) 



Here, $7,511 is the observed 511 keV gamma-ray fiux and 

Rqc is the distance of the Earth from the galactic center. 
Inserting the observed fiux (^>7,5ii — 10"'^ ph/cm^s) and 
Rqc ~ 2.5 X 10^^ cm, we arrive at 

Tdm 4 X 10^*^ seconds/mdm(MeV), (5) 

which is considerably longer than the age of the universe. 

III. AXINO DECAY 

Via R-parity violating couplings, axino LSPs may de- 
cay to Standard Model particles. The decay width for 
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such processes depends on the axion model and the na- 
ture of the R-parity violation considered. 

There are two classes of invisible axion models, KSVZ 
[29] and DFSZ [30]. Both of these introduce at least 
one extra field, which breaks the Peccei-Quinn (PQ) 
symmetry at a high scale, /„. Supersymmetric versions 
of these models give rise to two corresponding classes of 
axinos. 

In the first of these scenarios (KSVZ), the field $ cou- 
ples to a pair of heavy quark states via the superpotential 
coupling, W = X^QQ. In the second scenario (DFSZ), 
Standard Model fermions carry PQ charge. However, 
they do not have direct couplings to the PQ field. The 
PQ field, however, couples to the Higgs sector which con- 
tains two Higgs doublets. In both scenarios, the axino has 
a gaugino-gauge boson coupling proportional to the PQ 
anomaly [21] 



as 



(6) 



CaYY is a model dependent number of 0(1). 

We will be mainly interested in axino-fermion-sfermion 
couplings, in particular axino- Icpton-slcpton couplings, 
g-if, in both of those scenarios. In the KSVZ sce- 
nario, there is no direct coupling between the axino and 
Standard Model fields since they arc not charged under 
U{l)pQ. However, such an coupling could be induced at 
loop-level via a Bino-^y-lepton loop [22] as 
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On the other hand, in the DFSZ scenario, since the axino 
has a higgsino component, there is going to be a direct 
coupling 
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fa 
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where g is of the size of the gauge coupling suppressed 
by the mixing angle between higgsinos and gauginos. For 
our estimate, it is useful to pull out the common factor 
in g^ij (with the rough order of magnitude identification 
V ~ Ml in mind) and write g-ij = gv/ fa- From Eq. 7 and 
Eq. 8, we can roughly estimate g ^ 10~^ in the DFSZ 
case and g ~ 10"'' in the KSVZ case. 

The exact decay width depends on the details of the 
model, such as the Higgs potential and the spectrum and 
mixings of the superpartners. Rather than going into a 
detailed study, we give here an order-of-magnitude esti- 
mate. The axino life time is estimated to be 
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A is the R-parity violating leptonic trilinear coupling 
which appears in the superpotential as = Xij^LiLjE^. 
Positrons could be produced in the decays a — > z/^e+e" 
or a ^ v^e^e~ which result from the couplings A311 
and A211, respectively. A211 is constrained by charge 
current universality as A211 ^ 0.1(me„/200 GeV). A311 
is constrained by T[t ^ evv) IT{t ^ jivv) as A311 < 
0.12(me„/200GeV). See Ref. [17] for details. For axinos 
with a mass in the range of 1-300 MeV, for a wide range 
of couplings it is possible to obtain the desired lifetime 
found in Eq. 5. 

In addition to trilinear R-parity violating decays, bilin- 
ear couplings of the form fiiLiH2 could also contribute 
to the 511 keV gamma-ray production. The impact of 
this term on the axino lifetime was studied in Ref. [19]. 
This bilinear term induces mixing between higgsinos (and 
hence photinos) and neutrinos. Therefore, axinos can de- 
cay via a ^ ^+1/. Hi is constrained by the diffuse gamma 
ray background to be < keV [19]. In this case, the axino 
lifetime is about 10^^ sec which is potentially interesting. 
However, since this decay produces 7-1-1^ rather than 
positrons, this mode is useful for explaining the 511 keV 
line structure only if the mass of the axino is precisely 
1022 keV. Such a scenario is indeed fine-tuned, although 
it is technically natural. Another potentially important 
mode induced by the bilinear coupling is a ^ t+ -|- 7r~ . 
However, it is kinematically forbidden for the range of 
axino masses we are interested in. 



IV. DISCUSSION AND CONCLUSIONS 

In this letter, we have demonstrated that a light axino 
(1-300 MeV), in either the KSVZ or DFSZ axion mod- 
els, with trilinear R-parity violating couplings could be 
responsible for the 511 keV line emission observed from 
the galactic bulge. In this scenario, axinos constitute the 
major component of the cold dark matter and are present 
in the galactic halo with a cusped distribution (7 ^ 1.2). 

At this time, we can not exclude the possibility that 
poorly understood conventional astrophysics is responsi- 
ble for the observed positron production in the galactic 
bulge. To differentiate such a scenario from more ex- 
otic sources, such as light decaying particles, future tests 
must be made [31]. Additionally, if decaying axinos are 
the source of the 511 keV emission, signatures of super- 
symmetry, and R-parity violation will likely be observed 
at the LHC. 

As this letter was being completed, an article appeared 
which also discussed decaying particles as the source of 
the observed 511 keV emission [32]. They considered 
decaying sterile neutrinos with rather constrained mixing 
parameters. They also discussed decaying scalars with 
gravitational strength interactions. 

Acknowledgments: DH is supported by the Leverhulme 
Trust. LW is supported in part by U.S. Department of 



3 



Energy under grant DE-FG02-95ER40896. We would 
like to thank D. Chung and L. Roszkowski for useful dis- 
cussions. 



[1] J. Knodlseder, et.al, Accepted for publication in A&A, 
arXiv:astro-ph/0309442; P. Jean et al, arXiv:astro- 

ph/0309484. 

[2] R. Kinzer, P. Milne, D. Kurfess, M. Strickman, 
W. Johnson and W. Purcell, Astrophys. J. 559, 282 
(2001). 

[3] P. A. Milne, J. D. Kurfess, R. L. Kinzer and M. D. Leis- 

ing, New Astron. Rev. 46, 553 (2002) [arXiv:astro- 
pli/0110442], and references therein. 

[4] R. E. Lingenfelter and R. Ramaty, Positron-Electron 
Pairs in Astrophysics, eds. M. L. Burns, A. K. Harding 
and R. Ramaty, AIP Conference Proceedings, 267. 

[5] R. Ramaty, B. Kozlovsky and R. E. Lingenfelter, Astro- 
phys. J. Suppl. 40, 487 (1979). 

[6] B. Kozlovsky, R. E. Lingenfelter and R. Ramaty, Astro- 
phys. J. 316, 801 (1987). 

[7] P. A. Sturrock, Astrophys. J. 164, 529 (1971). 

[8] P. A. Milne, L. S. The and M. D. Leising, arXivrastro- 
ph/0104185. 

[9] M. Casse, B. Cordier, J. Paul and S. Schanne, 

arXiv:astro-ph/0309824. 
[10] S. E. Woosley and A. Heger, arXiv:astro-ph/0309165. 
[11] J. Bland-Hawthorn and M. Cohen, Astrophys. J. 582, 

246 (2003) [arXiv:astro-ph/0208553]. 
[12] M. Pohl, arXiv:astro-ph/9807268; C. D. Dermer and 

J. G. Skibo, arXiv:astro-ph/9705070. 
[13] C. Boehm, D. Hooper, J. Silk and M. Casse, arXiv:astro- 

ph/0309686. 

[14] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977). 

[15] P. Fayot, Phys. Lett. B 69, 489 (1977). 

[16] S. Weinberg, Phys. Rev. D 26, 287 (1982); N. Sakai and 

T. Yanagida, Nucl. Phys. B 197, 533 (1982). 
[17] H. K. Dreiner, arXiv:hep-ph/9707435. B. C. Allanach, 

A. Dedes and H. K. Dreiner, Phys. Rev. D 60, 075014 

(1999) [arXiv:hep-ph/9906209]. B. C. Allanach, A. Dedes 
and H. K. Dreiner, arXiv:hep-ph/0309196. 

[18] S. A. Bonometto, F. Gabbiani and A. Masiero, Phys. 
Rev. D 49, 3918 (1994) [arXiv:hep-ph/9305237]; T. Goto 
and M. Yamaguchi, Phys. Lett. B 276, 103 (1992); 
T. Asaka and T. Yanagida, Phys. Lett. B 494, 297 

(2000) [arXiv:hep-ph/0006211]; J. E. Kim, arXiv:astro- 
ph/9711310. K. Rajagopal, M. S. Turner and F. Wilczek, 
Nucl. Phys. B 358, 447 (1991); T. Goto and M. Yam- 
aguchi, Phys. Lett. B 276, 103 (1992). 

[19] H. B. Kim and J. E. Kim, Phys. Lett. B 527, 18 (2002) 

[arXiv:hep-ph/0108101]. 
[20] E. J. Chun and H. B. Kim, Phys. Rev. D 60, 095006 

(1999) [arXiv:hep-ph/9906392]. 
[21] L. Covi, J. E. Kim and L. Roszkowski, Phys. Rev. Lett. 

82, 4180 (1999) [arXiv:hep-ph/9905212]; 
[22] L. Covi, H. B. Kim, J. E. Kim and L. Roszkowski, JHEP 

0105, 033 (2001) [arXiv:hep-ph/0101009]; 



[23] L. Covi, L. Roszkowski and M. Small, JHEP 0207, 023 

(2002) [arXiv:hep-ph/0206119]. 

[24] T. P. Walker, G. Steigman, D. N. Schramm, K. A. Olive 
and H. S. Kang, Astrophys. J. 376, 51 (1991). 

[25] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. 
B 303, 289 (1993). 

[26] J. R. EUis, K. A. Olive, Y. Santoso and V. Spanos, 
arXiv:hcp-ph/0312262; J. L. Feng, A. Rajaraman 
and F. Takayama, Phys. Rev. D 68, 063504 (2003) 
[arXiv:hep-ph/0306024]; R. H. Cyburt, J. R. Ellis, 
B. D. Fields and K. A. Olive, Phys. Rev. D 67, 103521 

(2003) [arXiv:astro-ph/0211258]. 

[27] T. Gherghetta, G. F. Giudice and A. Riotto, Phys. Lett. 

B 446, 28 (1999) [arXiv:hep-ph/9808401]. 
[28] J. F. Navarro et al., arXiv:astro-ph/0311231; J. Die- 

mand, B. Moore and J. Stadel, arXiv:astro-ph/0402267. 
[29] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shif- 

man, A. L Vainshtein and V. L Zakharov, Nucl. Phys. B 

166, 493 (1980). 
[30] M. Dine, W. Fischlor and M. Srednicki, Phys. Lett. B 

104, 199 (1981); A. R. Zliitnitsky, Sov. J. Nucl. Phys. 

31 (1980) 260 [Yad. Fiz. 31 (1980) 497]. 
[31] D. Hooper, F. Ferrer, C. Boehm, J. Silk, J. Paul, 

N. W. Evans and M. Casse, arXiv:astro-ph/0311150. 
[32] C. Picciotto and M. Pospelov, arXiv:hep-ph/0402178. 



4 



